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ABSTRACT 

We present the first spectroscopic metalhcities of three M dwarfs with known 
or candidate planetary mass companions. We have analyzed high resolution, 
high signal-to-noise spectra of these stars which we obtained at McDonald Ob- 
servatory. Our analysis technique is based on spectral synthesis of atomic and 
molecular features using recently revised cool-star model atmospheres and spec- 
trum synthesis code. The technique has been shown to yield results consistent 
with the analyses of solar-type stars and allows measurements of M dwarf [M/H] 
values to 0.12 dex precision. From our analysis, we find [M/H] = -0.12, -0.32, and 
-0.33 for GJ 876, GJ 436, and GJ 581 respectively. These three M dwarf planet 
hosts have sub-solar metallicities, a surprising departure from the trend observed 
in FGK-type stars. This study is the first part of our ongoing work to determine 
the metallicities of the M dwarfs included in the McDonald Observatory planet 
search program. 

Subject headings: planetary systems - stars: abundances - stars: late-type - 
stars: individual (GJ 876, GJ 436, GJ 581) 



1. INTRODUCTION 



The detection of the fir st candidate planetary mass companion orbiting a solar-type 
star (IMayor fc Oueloz Ill995l ) ushered in the era of extrasolar planet research in astronomy. 
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Tremendous progress in this area has been made in the decade since and more than 200 can- 
didate planets have been announced. Despite the success of the current detection methods, 
unambiguous direct imaging of an extrasolar planet orbiting a star is still an elusive goal. 
Therefore, most of the knowledge regarding planet formation and evolution that has been 
garnered from the systems detected so far has come from statistical studies of the system 
properties and the host stars themselves. One of the interesting findings from these studies 
is the trend toward higher photospheric metal abundances in extrasolar planet hosts stars 
relative to stars without detected planets. 



Gonzalez I (119971 ) first noted the high metallicities of the first four stars which were 
found to exhibit radial velocity variations attributable to a planetary mass companion. This 
trend was found to continue as more stars were identified as potential extrasolar planet 



hosts with the Doppler method, and fol 



(Fuhrmann et al. 


1997 




1998; 


Gonzalez 


2000. 


2003. 


2004 


) - 


2005; 


Laws et al. 


2003 



owed up wi t h high precision abun d ance analyses 



1998 



1999 



Gonza. 



Fischer fc Valenti 



ez et al. 2001; Santos et al. 



20051 ). The most likely expla- 



nation for this observed trend is the so-called "primordial" hypothesis. That is, the high 
photospheric metal abundances in the h ost stars are relics of protostellar clouds and disks 
with a proportionally high metal content (ISantos et al. II2004J : iFischer fc Valenti II2005I ). It is 
theorized that high-mass planet formati on is increased in hig h metal-content protoplanetary 
disks under the core-accretion paradigm (jPoUack et al. Ill996l ). This hypothesis explains why 
more Jupiter and higher-mass planets have been detected around stars with high metallici- 
ties. 

The majority of extrasolar planets have been found around FGK-type stars as these are 
the stellar types that make up the majority of targets in Doppler surveys. High precision 
abundance analyses for these types of stars are relatively straightforward and, therefore, 
these are the types of stars for which the metallicity - giant planet connection h as been 



established. However, the majority of stars in the solar neighborhood are M dwarfs (IHenry 



19981 ) and a complete understanding of planet formation must necessarily include late-type 
stars. 

To date, only one M dwarf, GJ 876, is known to harb or a Jupiter-mass companion 



(IDelfosse et al.lll998l ; iMarcy et al. Ill998l ; [Benedict et al.ll2002l ). In addition to the astromet- 



rically confirmed outer planet, another Jupiter- mass pla net and a very low-mass planet in 



shorter-period orbits have been detected around GJ 876 (IMarcy et al. Il200ll : iRivera et al. 



2005 



making it also th e only M dwa rf multi-planet host known to date. The M dwarfs GJ 



436 (IButler et al.ll2004j ) and GJ 581 (IBonfils et al.ll2005bf ) are hosts to candidate Neptune 



mass planets in short-period orbits based on their radial velocity variations. Other planets, 
including a gas giant planet and a ~ 5.5 A^© planet, have been detected around suspected M 
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dwarfs using the microlensing technique (iBeauheu et ah I l2006l : iGould et ah 1 12006| ) . These 
host stars are still confused with the source stars and unavailable for further study. The first 
star that was identified as an extrasolar plaii et host with the microlensing technique was 
thought to be an M d warf (IBond et al. 1120041) . but recent observations have shown that it 



is actually a K dwarf (IBennett et al 



20061 



In this Letter, we present the results of our abundance analysis of three M dwarfs with 
a known (GJ 876) or candidate (GJ 436 and GJ 581) planetary mass companion. In §2 we 
describe our high resolution (R) and signal-to- noise (S/N) spectroscopic observations. We 
discuss our analysis and present our results in §3. We briefly discuss our results and ongoing 
effort to determine the metallicities for all the M dwarfs being monitored in a Doppler survey 
4. 



m 



OBSERVATIONS AND DATA REDUCTION 



We observed GJ 876 and 436 using the 2.7 m Harlan J. Smith telescope at McDonald 
Observatory on November 20, 2003 and January 24, 2005 respectively. Data were obtained 
with the 2dcoude spectrograph (ITuU et al. Ill995l ) equipped with a 79 gr mm~^ echelle grating 
and 8.2"x 1.2" slit. Two 30 minute exposures were taken for both objects and co-added before 
order extraction. 

GJ 581 was observed using the 9.2 m effective aperture Hobby-Eberly Telescope (HET) 
at McDona ld Observatory on May 11, 2006 with the High Resolution Spectrograph (HRS 
Tulllll998h fed by a 2"optical fiber. The HRS was used in the "R = 60,000" mode with a 
316 gr mm~^ cross-dispersion grating. The cross- dispersion grating was positioned so that 
the break between the two CCD chips was at 7940 A. Two 10 minute exposures were taken 
for GJ 581 and co-added before order extraction. 



CCD reduction and optini al order extraction were carried out using the REDUCE pack- 
age (jPiskunov fc Valenti 1120021 ). The wavelength calibrations for each object were calculated 
based on the identification of roughly 1000 lines in thorium-argon emission spectra taken 
at the beginning of each respective night and have RMS precisions of 0.002 A. The final 
one-dimensional spectra of GJ 876, 436, and 581 have S/N, of 430, 360, and 190 pixel"^ 
respectively at 8700 A. The measured resolving powers were roughly 50,000 for the 2.7m and 
60,000 for the HET data. 
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ANALYSIS AND RESULTS 



We ana lyzed the ob s erved spectra of the M dwarf planet hosts using the t echnique 
described by lBean et al. I (120061 ) to determine their metallicities. iBean et al. I (120061 ) utilized 
analyses of both components of solar-type and M dwarf visual binaries to test and improve 
their spectroscopic analysis technique and cool-star model atmospheres. Their test was 
based on the assumption t hat un evolved stars in bound systems have the same photospheric 
abundances. iBean et al. I (120061 ) showed that their method yields metallicities for M dwarfs 
consistent with the results given by standard analysis techniques applied to solar-type stars. 
We briefly describe the technique here, and refer the reader to that paper for a complete 
description. 

Our analysis of the M dwarf planet hosts relied on fitting synthetic spectra to their 
observed spectra. We used spectral regions containing a strong TiO bandhead and rela- 
tively clean atomic line profiles as the constraints and the statistic as our goodness of 
fit metric. We generated synthetic spectra with an updated version o f the plane-pa rallel, 
local thermodynamic equilibrium (LTE), stellar analysis code MOOG (jSneden 1119731 ). We 
adopted a grid of model atmospheres which were cor nputed for this particul ar purpose using 
the model atmosphere code PHOENIX (version 13. iHauschildt et al.lll999l ). We generated 
a model atmosphere with arbitrary parameters to be input into MOOG by interpolating in 
this grid. 

We determined the stellar parameters effective temperature, Tg//, metallicity, [M/H0, 
and microturbulent, ^, and macroturbulent, 7], velocities directly from the spectral fitting. 
We constrained the stellar surface gravity using an empirical log g - Ai relationship that 
was derived from recent measurements of M dwarf radii. We estimated the masses of the 
objects to use with this relationship based on the - M. relationship in iDelfosse et al. 



(120001 ) and using the K magnit udes and parallax es from the Two Micron All Sky Survey 



(2MASS) point sourc e catalog (ICutri et al.l 120031 ) and the Hipparcos catalo g (lESA 



Ferryman et al.l 119971 ) respectively. We assumed the abundances given by lAsplund et al. 



1997 



(120051 ) as the reference solar abundances and linear alpha element (O, Mg, Si, Ca, and 
Ti) and C enhancement r elationships as functions of [Fe/H] from the data presented by 
AUende Frieto et al.l J2004h . 



The results from our analysis of the M dwarf planet hosting stars GJ 876, 436, and 581 
are given in Table [H We find [M/H] = -0.12, -0.32, and -0.33 for GJ 876, GJ 436, and 



^We adopt the standard spectroscopic notation: for elements X and Y, log e(X) = logiQ{Nx /Nh) + 12.0, 
[X/Y] = \ogio{Nx /Ny)* - \ogio{Nx /Ny)q, and Nx is the number density of element X. 
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G J 581 respectively. We adopt the standard uncertainties in the parameters derived using 



this technique as given by iBean et al. I (120061 ) . They are 48 K, 0.10 dex, 0.12 dex, 0.15 km 



s ^, and 0.20 km s ^ for Te//, log g, [M/H], ^, and r] respectively. These uncertainties 
which we consider to be inherent in the analysis technique, were calculated based on the 



agreem ent between the M dwarf and solar-type star visual binary pairs in the iBean et al. 



( I2OO6I ) analysis. This external error estimation method includes correlated uncertainties from 
all the other determined parameters and yields uncertainties that are more conservative and 
robust than errors derived from only considering the uncertainty in the fitting process from 
the spectral S/N. 

Plots of the observed spectrum and best fit synthetic spectrum for GJ 876 are shown 
in Figures [T] and [21 Synthetic spectra computed with [M/H] values 0.3 dex lower and 
higher than the best fit value are also included in these Figures to illustrate the sensitivity 
of our measurement technique. The "high" and "low" metallicity synthetic spectra clearly 
do not match the observed spectrum as well as the synthetic spectrum computed with the 
determined stellar parameters. 

An interesting aspect of this result is the closeness of the derived Teff values for all three 
stars (range of 20 K) despite a range of 1.5 spectral types. The explanation for this is that 
metallicity and effective temperature are degenerate in the M dwarf spectral classification 
system. Therefore, only a detailed analysis such as the one we have employed can break the 
degeneracy and give a precise estimate of these parameters for an M dwarf. 



DISCUSSION 



Butler et al.l (120041 ) report an occurrence rate of Jupiter-mass planets (0.5 A4 jup < 
Ai <13 J^ .hi,p) with orbital se mi-major axes, a, < 1 AU of 3.5% around FGK-type stars. 



In contrast, lEndl et al. I (120061 ) found a frequency of 0.46% with an upper limit of 1.27% for 



the occurrence of Jupiter-mass planets around M dwarfs ba sed on a dedicated survey of the 



spectral type. This result is similar to that also reported by [But 



■er et al.l (12004 0.7 %) based 



on a survey of 150 M dwarfs with potential overlap with the lEndl et al. I (120061 ) sample. 
There have been no detections of so-called "hot Jupiters" (a ~ 0.04 AU) around M dwarfs 
despite the stronger sensitivity of the Doppler de tection method to t hese types of planets 
around low-mass stars. Around FGK-type stars, iMarcy et al. I (120051 ) cites a frequency of 
1.2% for hot Jupiters. 



While the limits that can be currently placed on the frequency of Jupiter-mass planets 
around M dwarfs are not entirely inconsistent with those of FGK-type stars, there does 
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seem to be a trend to fewer high-mass planet detections around M dwarfs. If further results 
support this observation, i t would be coii s istent with the predict i ons of the core accretion 
planet formation model. iLaughlin et al.l (120041 ) and llda fc Lin I (120051 ) have shown that 
the formati on pro b ability of high-mass planets decreases with stellar mass with this model. 
Conversely, IBoss I (120061 ) suggests that giant planets might actually form more efficiently 
around M dwarfs if the gravitational instability mechanism is considered. 

Further clouding the issue is the question of host star metallicity. As mentioned earlier, 
planets are more often detected around stars with high metallicities. In this Letter we have 
presented the results from our spectroscopic metallicity analysis of three M dwarfs that 
harbor extrasolar planets. We find that all three have sub-solar metallicities which is a 
departur e from the observed t rend in the FGK-type stars that harbor extrasolar planets. In 
contrast. iBonfils et al.l (l2005al ) presented metallicity measurements for these M dwarfs based 



on a lower precision photometric relationship. Our derived [M/H] values are lower by 0.09, 
0.08, and 0.34 dex fo r GJ 876, GJ 581, and GJ 436 respectively than those determined by 
Bonfils et al.l (j2005al ). In the case of the first two, the values are well within the overlapping 
errors for the two measurements (0.32 dex), while our measurement for GJ 436 is just outside 
this differential range. 



Taken together, the results from our analysis and that of lBonfils et al.l (l2005al ) do appear 
to rule out super-solar metallicities for these M dwarf planet hosts. This result raises some 
interesting questions. Are the metallicities for these stars representative of the metallicities 
of the M dwarfs on planet search programs and might that explain the lower detection rates 
of planets for the M dwarfs? If that were the case, are the solar neighborhood M dwarfs in 
general metal deficient relative to the other spectral type? Or, what is causing the selection 
effect to lower metallicity M dwarfs for the planet search programs? 

We plan to follow up this preliminary result by applying the same spectroscopic metal- 
licity analysi s technique to the M dwarfs that have been surveyed for extrasolar planets as 



described by lEndl et al. I (120061 ). The results from a larger sample should make it easier to 



disentangle the effects of stellar mass and metallicity on planet formation. 
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Table 1. Spectral types and derived stellar parameters for the M dwarf planet hosts. 



Name 


Spectral Type 




log g 


[M/H] 




■n 






(K) 


(cgs) 




(km s-l) 


(km s-1) 


GJ 876 


M4 


3478 


4.89 


-0.12 


0.77 


0.64 


GJ 436 


M2.5 


3498 


4.80 


-0.32 


1.02 


0.00 


GJ 581 


M3 


3480 


4.92 


-0.33 


0.91 


1.35 



''Adopted uncertainties are 48 K, 0.10 dex, 0.12 dex, 0.15 km s~^, and 0.20 
km for the derived parameters T^yy, log g, [M/H], ^, and r] respectively. 
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Fig. 1. — Spectral region near the strong TiO 7 R2 - bandhead for GJ 876 (histogram). 
The best fit used to determine the stellar parameters is over-plotted (solid black line). For 
comparison, synthetic spectra computed with [M/H] values 0.3 dex lower (dotted line, solid 
red in the electronic edition) and higher (dashed line, solid blue in the electronic edition) 
than the best fit value are also over-plotted. 
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Fig. 2. — Fit of synthetic spectra (solid line) to atomic line profiles (points) for GJ 876. The 
filled points were used in the fitting process; the open points were ignored. For comparison, 
synthetic spectra computed with [M/H] values 0.3 dex lower (dotted fine, solid red in the 
electronic edition) and higher (dashed line, solid blue in the electronic edition) than the best 
fit value are also over-plotted. The panels are sorted by wavelength and the linear scaling in 
both parameters is the same throughout. The lines in each half, top and bottom, make up 
a contiguous spectral order in our observed spectra. All apparent "lines" in the figure that 
aren't fit are actually multiple TiO lines. 



